Subarea law of entanglement in nodal fermionic systems 
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We investigate the subarea law scaling behavior of the block entropy in bipartite fermionic systems 
which do not have a finite Fermi surface. It is found that in gapped regimes the leading subarea 
term is a negative constant, whereas in critical regimes with point nodes the leading subarea law 
is a logarithmic additive term. At the phase boundary that separates the critical and non-critical 
regimes, the subarea scaling shows power-law behavior. 
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The study of entanglement across quantum phase tran- 
sitions has attracted substantial activity in the physics 
community. This interest was initially boosted by the 
discovery of a connection between the scaling properties 
of the concurrence, which is a measure of entanglement, 
and the quantum phase transition in the one-dimensional 
(ID) transverse Ising model[l„ J,]. Subsequently, the scal- 
ing properties of the von Neumann entanglement entropy 
(EE), which is defined as Sl = — Trpi^d log pi^d of its 
reduced density matrix p^d = Trif\i^d\^!){^\ and L is 
the length of a contiguous system block L'^, show a dra- 
matic difference between ID critical and non-critical sys- 
tems: in critical regimes it exhibits logarithmic scaling 
Sl oc logL, whereas in non-critical regimes it saturates 
to a constant as L ^ oo Conformal field the- 

ory (GET) tells us that the analytical form of the EE 
in ID critical cases is given hy Sl = log2 L, where 
c and c are the holomorphic and antiholomorphic cen- 
tral charges which are universal factors for a given class 
of critical systems 5, 6]. In this sense, if the correlation 
length diverges, the EE exhibits logarithmic scaling with 
block size, and if the correlation length is finite, we sim- 
ply obtain a constant EE. 

In higher dimensions, d > 1, the scaling behavior of 
the EE is complicated and far from fully settled. On 
one hand, in bosonic systems it has been shown that the 
area law Sl ~ L'''~^ hold for both critical and non-critical 
cases. On the other hand, in fermionic systems, the scal- 
ing behavior of EE highly depends on the topology of 
the Fermi surface: the area law was shown to be violated 
in systems possess a finite Fermi surface 0, S]; but it 
still holds in gapped systems and systems having point 
nodes [13, [ill. Il2|. It is remarkable that in systems 
where the area law holds, the link between correlation 
length and the scaling behavior of the EE is broken: the 
leading term of EE fails to tell if the correlation length 
diverges in the system. This is not too surprising when 
the leading term follows the area law since an area law 
only reflects short-range correlations in the system. How- 
ever, this motivates us to study the subarea behavior of 
the EE, where the link between correlation length and 
EE may still be active. 



The subarea behavior of EE has been investigated in 
some 2D systems. Noticeably, it exhibits some very in- 
teresting properties. In non-critical systems, the com- 
mon expectation for the subleading behavior is a con- 
stant, which in some cases can be used to characterize 
topological order 1^ [13, [i3| • For critical systems in 
the universality class of z=2 conformal quantum critical 
points, there exists a universal logarithmically divergent 
correction to the non-universal area law, which is deter- 
mined by the shape of the partition and by the central 
charge. Very recently, a universal logarithmically diver- 
gent correction in the 2D random transverse Ising model 
was also observed [3]. Following their approach, one 
might expect that a general logarithmic additive correc- 
tion generally exists in critical systems. 

In this paper, we examine the behavior of the sublead- 
ing term of the EE in a 2D quadratic fermionic Hamil- 
tonian. Our results show that in the non-critical regimes 
the leading subarea term is a negative constant, propor- 
tional to the square root of the correlation length. On 
the other hand, in quantum critical phases which have 
only point nodes in their excitation spectrum the leading 
subarea law is found to be a logarithmic additive term. It 
is interesting to note that the subarea law shows unusual 
behavior in the vicinity of the phase boundary separat- 
ing the critical and noncritical regimes. According to the 
careful numerical analysis discussed below, this unusual 
behavior follows a nonadditive power law with a < 1. 

The particular Hamiltonian we consider in this paper 
is a bilinear spinless fermionic model on a 2D lattice 
with pairing interaction between nearest-neighbor lattice 
sites [lO|, 



H = 



(ij) 



(1) 



Here A is the chemical potential, and 7 is a p-wave pair- 
ing interaction. The sum J2{ij) extends over nearest- 
neighbor pairs. Depending on the parameters 7 and A, 
this system has a rich phase diagram, including metallic, 
insulating and p-wave superconducting regimes. In previ- 
ous work, we distinguished three different phases whose 
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FIG. 1: Scaling behavior of the entanglement entropy Sl for a 
cut through the phase diagram at fixed interaction potential 
7 = 1.0. The chemical potential is varied in the window 
A £ [0.5,3.0], thus crossing through phases II and III which 
do not have finite Fermi surfaces. The parameters shown in 
the legend correspond to (A, 7). The inset shows the phase 
diagram of the 2D spinless fermion model (Eqn. 1) . 



FIG. 2: (a) Scaling behavior of S'f"" in phase II. (b) Scaling 
behavior of S'f in phase III. (c) The relation between the 



constant term const and ■ 
tential 7 — 0.5,7 = 1-0,7 



= for fixed finite interaction po- 
2.0 in phase III. (d) Dependence 

1 



of ^ v.s. L(A — 2) with the fixed interaction potential 7 
in phase III. 



different signatures are reflected by the scaling proper- 
ties of the leading term of the EE. Phase I is the case 
7 = and < A < 2. A finite Fermi surface (i.e. line 
nodes) exists in this regime. Phase II is the case with 
{0 < A < 2,7 > 0}. In this regime, only point nodes 
exist in the excitation spectrum. Phase III {A > 2} is an 
insulating state characterized by a gap in the excitation 
spectrum. The resulting phase diagram is shown in the 
inset of Fig. [1] The ground state in this model is known 
to be a Gaussian state whose reduced density matrix can 
be obtained exactly with the usa ge o f Grassman algebra 
and a Bogoliubov transformation^^ [2lj . Therefore, one 
can calculate the block entropy Sl directly, in the ther- 
modynamic limit, using the correlation matrix[ia El. 
Previously fio'l, it was found that in phase I the leading 
scaling behavior of the EE is Sl ^ LlogL and the coef- 
ficients of this term are well described by the analytical 
form deduced by the Widom conjecture[8|. In phase III, 
the area law Sl ^ L holds. Interestingly, in phase II 
area law behavior is also observed (Fig. [T|) although the 
system has a diverging correlation length in this case. 

To study the subarea law of the EE, we need to first 
subtract the dominant area law contribution. Since we 
are able to compute the block entropy Sl directly in 
the thermodynamic limit, the dominant area law con- 
tribution depends linearly on the size L of the block: 
Sl ~ clL, for a two-dimensional system. This simple 
linear dependence should be contrasted with the domi- 
nant, non-monotonic dependence of EE on block size in 



a finite-size system 15|, Il6l . [17[. For example, the dom- 
inant contribution to the EE of a block of size L in a 



one-dimensional system with Hamiltonian (JlJ has the fol- 
lowing non-monotic form Q: 



5.^-log(-sm-). 



(2) 



00, this takes the 
Similarly, we find a 



In the thermodynamic limit N 
monotonic form 5*^ = 1 logj L + . 
monotonic, dominant contribution to the block entropy 
in two dimensions, which allows us to extract sub-leading 
terms reliably. 

We will focus on the quantity S"! = LSl+i - (i + 
1)Sl, and we will also examine the second-order deriva- 
tive ^1^- It is then straightforward to see that if 
the EE follows a behavior Sl ^ aL — const + 0(-^), 

sl"'' const and ^ Oi^w)- In contrast, if the 

EE scales as Sl ^ aL - blogL + 0(1), then S'f"'' ~ 



-5Llog(L + 1) + 5(L + l)logL - 61ogL and ^ - 
for L ^ 1. Interestingly, if the EE has a power-law sub- 
leading correction, 5*^ = aL — bL" + 0(1) with a < 1, 
5f"'' ~ 6L" and ^ - j^. 

A significant difference of S'f between phase II and 
phase III is observed despite the fact that both phases 
obey an area law in leading order. The results are shown 

in Fig. H We find that in phase II i'f blogL, with 

6 > (Figl2{a)). In contrast, in phase III, one observes 
that S'f quickly converges to a positive constant as L 
increases (FiglJJb)). These results are consistent with 
the following scaling forms for the entanglement entropy: 
Sl — aL—blogL+. . ., with 6 > in phase II, and in phase 
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III Sl = aL — const + . . ., with const is an positive non- 
universal constant which depends on (A, 7) . As noted 



in Ref. [ij], the sign of a sub- leading correction to the 



entropy may be positive or negative, since the leading 
term insures that the entropy is positive, as it should 
be. In fact, both the conformal quantum critical points 
studied in Ref. 14 1 and the critical points in Ref. (loj 



have subleading logarithmic corrections with the same 
sign as found in phase II. 

In contrast to the fixed points described in Refs. 



14 . Il9l | , the z = 1 critical phase II has a non- universal 
sub-leading correction [2^. For example, in Figl2ja), we 
see that moving around within phase II changes the slope 
of S'f"'' dramatically. Similarly, the sub-leading correc- 
tion is non-universal in phase III, as seen in Figl^^b). 
This is a non-surprising result, considering that the phase 
is gapped with a finite correlation length ^. Interestingly, 
the absolute value of this correction follows the same 
trend as the correlation length ^ in phase III, namely 
it increases with increasing pairing interaction 7 and de- 
creases if A increases. Therefore, let us now focus on the 
relationship between the constant term of S'f"'' in phase 
III and the parameters (A, 7). In Figl2^c), we examine 
the relation between the constant term and ^-t fixed 

values of the interaction potential 7 = 0.5, 1.0, 2.0. One 
finds good linear fits and observes that the slopes of these 
lines increase with 7, which suggests the following form: 
SL^aL- in phase III. FigE^d) gives the relation 

between ^ and L{X — 2). For fixed 7, one finds that 
all curves collapse. With the the band gap Aq = A — 2 
and the knowledge that the transition from phase II to 
phase III has z = 1 criticality, it is plausible to infer 
from the above graph that the subarea law in this regime 
has a uniform finite-size scaling form: Sf"** ~ fi^j)^ 

with f{x) = when x — > 00, and f{x) — const, when 
a; ^ 0. 

The critical scaling regime of the boundary between 
phase II and phase III provides a striking, universal cor- 
rection to the area law (Figl3]). Here the dominant sub- 
leading behavior is not logarithmic (See Fig[3l^a)). We 
find that S'f exhibits significant curvature as a function 
of log L in phase III, in contrast to the linear behavior 
in phase II (denoted by filled circles in FiglSfa)). Fur- 
thermore, using the log-log graph (Fig[3Jb)), one finds 
that S'f follows a power-law relation over almost two 
decades of scaling. Finally, the exponent we extract from 
fitting is close to a = 0.5 (eg: a = 0.493 ± 0.001 for 
A = 2.0 ,7 = 2.0). In fact, all the curves can be fit with 
the log-linear law logy = .blogx + b(pj), with a slope that 
is independent of 7, as long as one remains on the critical 
line separating phase II and phase III. This insensitivity 
of the power-law exponent is consistent with universal 
behavior of this sub-leading contribution to the entropy. 
By dimensional analysis, the coefficient b depends on the 
lattice spacing, which is consistent with the observed 7 
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FIG. 3: (a) Semi-log (Y/logX) graph and (b) Log-log graph 
of the scaling behavior of Sf in phase II and its bound- 
ary points, (c) Scaling behavior of in phase II and its 
boundary points. 



dependence of h in our numerical results. We believe this 
is the first example of a power-law correction to an area 
law observed in any dimension. 

As we mentioned above, the second-order differentia- 
tion 



of Sl can also be used to detect the scaling be- 
havior of the subarea law contributions to ■ This anal- 
ysis is performed in Fig. EI^c)). One finds that fol- 
lows a perfect inverse power law relation. However, there 
exist obvious differences between phase II and the bound- 
ary region, where the exponent of this inverse power law 
is not the same. For phase II, after fitting the curves in a 
log-log graph, we obtain (3 very close to —2.0. A typical 
example with A = 1.0, 7 = 1.0 is shown in Fig. [2Ic), the 
simple linear fitting gives us /3 = 1.988±0.001 in this case. 
This result suggests Sl ^ aL — h log L and is consistent 
with the results from the above analysis of in Fig. 

In the boundary regime, one finds the power factor 
of ^ to be very close to -1.5(eg: (3 = -1.509 ± 0.002 
for A — 2.0 ,7 = 2.0), consistent with the above obser- 
vation S'^"^ ^ L^^-^ for all the four cases we studied. 
Therefore, both two numerical methods suggest unusual 
power law corrections in the subarea term of the EE. 

It is worthy to consider alternative possibilities for the 
subarea law in the boundary region apart from power-law 
corrections. One possibility is a log log L correction as a 
subleading term. In this case, would decay more 

rapidly than L~^, and so it cannot explain the results 
in Fig. [3fc) since in this graph is well described 

by [a ^ 0.5) in the boundary regime. A second 

possibility would be a term scaling as — log^ L. In this 



situation, it is true that 



dl? 



1 
I? 



can give 
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TABLE I: Coefficients of the Sa ■ Sa ^ aA ~ blnA + const 
in phase II (a) A = 1.0, 7 = 1.0. (b) A = 1.0, 7 = 2.0. 



Shape 



const 



Square (a) 0.224170(4) 0.1189(4) -0.191(1) 

Isosceles Righttriangle (a) 0.253410(6) 0.2431(8) -0.003(3) 

Square (b) 0.26230(5) 0.120(3) -0.33(1) 

Isosceles Righttriangle (b) 0.28719(3) 0.199(3) -0.201(9) 



us illusion of power law relation L~^+" when L is not 
too large. However, one finds that S'f"^ does not show 
a good linear relation on a log-log scale, and therefore is 
not consistent with the results in Fig. ^h). In brief, the 
unusual non-additive type subarea law is the most 
plausible leading candidate according to the above anal- 
ysis. Of course, the most straightforward evidence comes 
from our observation S'f 



L^fij) (fix) = const, at 
X 0) in the gapped phase, as shown in Figl2l[d). If 
one approaches the phase boundary, ^ — > 00, the form of 
subarea law should be S'f ~ L2 in analogy to the ID 
case where the finite size scaling form Sl ~ log'C in the 
non-critical phase changes to Sl ^ logi in the critical 
phase in the limit L » 1. It is noted that one cannot 
exclude the possibility that log log i , log^ L, even logi 
as the second leading order subarea-law term. 

The shape dependence of the logarithmic correction is 
another way to distinguish the origin of the non-universal 
corrections in phase II from the universal correction at 
the conformal quantum critical points IJ] . Here we com- 
pare results between square partitions and isosceles right- 
triangle (IR) partitions, the latter is taken as half of a 
square partition in our numerical simulation. The re- 
sults are presented in FigH) We find that the choice of 
shapes does not change the scaling behavior of 5f i.e. 
for phase II the subleading term is still logarithmic diver- 
gent (Figll^a)) despite that the prefactor of this logarith- 
mic term increases significantly in the IR case, and in the 
boundary region between phases II and III S'f still fol- 
lows power- law scaling behavior (FigdJ^b)). Currently, it 
is not clear whether under the change of geometry shape 
the power factor will change or not because the second- 
leading subarea law of Sl is not known and the inclusion 
of log log L or log L type terms can make the fits change 
significantly. 

In fact, one can fit our data using 5*^ =, Sa ~ o,A — 
blnA+c to obtain the coefficient directly for phase II. It is 
noted that the boundary length A for a square partition 
is 4L and for an IR partition, A = {2 + V2)L. For both 
of these two cases studied in Table I, it is shown that all 
coefficients are changing with the partition shape, which 
suggests that there does not exist simple relations for 
different shapes unlike the case in the 2D z — 2 conformal 
quantum critical points [l3|. Interestingly, for the case 
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FIG. 4: Comparison of Sl^ for two different partitions, (a) 
The phase II : A = 1.0, 7 — 1.0. (b) The boundary regime : 
A = 2.0, 7 = 1.0 



(A — 1.0, 7 — 1.0), we find that the prefactor b is equal 
to 0.1189(4) if the partition is a square and the prefactor 
b is equal to 0.2431(8) if the partition is an IR. This result 
is very similar to the prediction in Ref. IJ] , as the subarea 
law ~ § ~ 0.111c and ~ 0.242c if we 

consider c = 1 here. However, this is only a special case 
since one can immediately see there exist large differences 



in the ratio of 



in other cases, for example, A = 1.0, 



7 = 2.0 shown in Table I. 

In summary, our results suggest that there exist sig- 
nificant differences in the leading subarea term of the 
EE between non-critical models and critical models with 
point nodes in 2D bipartite fermionic systems. In the 
non-critical phases the leading order correction is a non- 
universal negative constant, proportional to the square 
root of correlation length. And in the critical phases with 
point nodes we find a subdominant logarithmic additive 
term with a non- universal coefficient. In the boundary re- 
gion between the point-node and the completely gapped 
phases the subarea law can not be described by a sim- 
ple logarithmic behavior anymore. Our analysis suggests 
that there is an unusual non-additive power law relation, 
which also depends on the shape of geometrical partition. 
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